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Abstract 



We revisit the results of recent electroweak baryogenesis calculations and 
include all allowed large CP-violating supersymmetric phases. If the phases 
are large, the resulting baryon asymmetry can be considerably larger than 
the observed value ub/s ~ 4 x 10^^"*^. Much of the asymmetry must therefore 
be washed out, and we argue that the upper bound on the light Higgs mass 
is larger than the value reported in previous work. 
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Typeset using REVT^ 



I. INTRODUCTION 



The universe is not matter-antimatter symmetric; as shown by primordial nucleosynthe- 
sis measurements, the baryon density to entropy density ratio ns/s is constrained to be 
about 4 X 10^^^ [|l|. The necessary ingredients for a theory to explain this asymmetry are 
the Sakharov criteria: baryon number violation, C and CP violation and non-equilibrium 
conditions . Electroweak baryogenesis is an attractive mechanism with the best chance of 
verification through traditional collider experiments 0. 

The central feature of electroweak baryogenesis is spontaneous symmetry breaking. At 
high temperatures, the vaccum expectation value of the Higgs field is zero but as the universe 
cools down over time, a second minimum appears in the potential at w 7^ 0. The universe 
undergoes a phase transition, from a "symmetric" phase where f = to the "broken" 
phase where f 7^ 0. If the phase transition is second order, which means that no potential 
barrier exists, then the universe is approximately in equilibrium throughout the transition 
and no baryogenesis can result (since the third Sakharov criterion is violated). If the phase 
transition is first order, then the phase transition occurs through quantum tunneling. As 
the universe cools, the probability of making a transition from symmetric to broken phase 
grows and when the transition occurs at a certain point in space, a bubble of broken phase 
forms and expands. Because of its first order nature, the phase transition does not occur 
throughout the whole universe at the same time, resulting in non-equilibrium conditions. 

The baryon asymmetry is generated as the wall of the expanding bubble passes through 
points in space. Particles in the unbroken phase close to the wall interact with the changing 
Higgs field profile in the wall and the presence of CP-violating couplings produces source 
terms for participating particles. Different chiralities couple with different strength when 
CP is violated and a difference occurs in the reflection and transmission probabilities for 
the two different chiralities. Due to rapid gauge, Yukawa and strong sphaleron interactions, 
the CP-violating source terms are translated into a net left handed weak doublet quark 
density which is finally converted into a baryon asymmetry by weak sphaleron decays. The 
asymmetry then diffuses through the bubble wall into the broken phase where the weak 
sphaleron interactions are exponentially suppressed. Subsequent washout of the baryon 
asymmetry can therefore be kept under control provided the first order phase transition is 
strong enough. 

The Standard Model satisfies the three Sakharov criteria, but it can not generate the 
required value of ub/s. This is because the only source of CP violation in the SM comes 
from the phase 5 in the CKM quark mixing matrix. The baryon asymmetry is ~ sin 5, but 
this value is suppressed by flavor mixing factors and the resulting asymmetry cannot be 
too large. Therefore, the phase transition has to be strongly first order to avoid washing out 
the produced baryon number, which translates into the criterion v{Tc)/Tc ^10. However, 
in the Standard Model, the phase transition is too weakly first order to prevent washout, 
unless the SM Higgs mass is less than 50 GeV, far below the present experimental limit [0 . 

The situation can be improved in the Minimal Supersymmetric Standard Model (MSSM) 
as the light right handed top squark contribution to the temperature dependent effective 
potential can push the value of v{Tc)/Tc up to acceptable values even for light Higgs masses 
allowed by experimental searches [p|- [T^JT6| -pO| . Also, supersymmetric extensions of the SM 
include additional sources, which in general can further enhance the possibilty of baryon 
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asymmetry production during the electroweak phase transition. As a result, a specific region 
in the MSSM parameter space corresponding to a heavy CP-odd Higgs boson, a hght CP- 
even Higgs boson, a hght right-handed stop and a heavy left-handed stop can provide a 
plausible framework for electroweak baryogenesis. 

The CP-violating interactions in the MSSM arise in the complex phases of the soft 
supersymmetry breaking terms in the Lagrangian and in the phase of fi. Since the largest 
contributions come from charginos and neutralinos [PTLPHl ^'^'^ ^^^^ importantly from the 
right handed stops, only the gaugino mass phases (pi, and (p^, the phase of and the 
phase of the stop trilinear parameter At are relevant for baryogenesis. These phases have 
to be small, typically < 10^^, if they are considered individually with sparticle masses 
0{TeV), otherwise they induce contributions to the electric dipole moments (EDMs) of the 
neutron and electron exceeding experimental limits pl|-p3|. However, these constraints can 
be avoided if relations among soft breaking parameters ensure cancellations of individual 
contributions to the EDMs P^ , P3| . In the most general case this possibility leads to models 
with light superpartner spectra and CP-violating phases of 0{1) ||3^. Remarkably, string 



motivated scenarios with non-universal gaugino masses can be described where such EDM 
cancellations occur naturally as a result of the SM embedding on five-branes . 



Different methods of calculating the baryon asymmetry, using both classical and quantum 
Boltzmann equations, have produced the required asymmetry, provided that sin*/?^ ~ 10~^ — 
10~^ p7|p0| -[3^ in agreement with the experimental limits on the EDMs of the electron and 



neutron, when the other MSSM phases are taken to be zero. In this paper we assume that 
large phases are allowed by the cancellation mechanism and at the same time the right 
handed stop is light enough to significantly modify the finite temperature Higgs potential. 
All other sfermions are heavier and their contribution to the Higgs potential is Boltzmann 
suppressed. 

We briefly review the generation of the baryon asymmetry and the phase dependence 
factorization emphasizing the significance of the + (fn phase combination. In the large 
phase scenario it is then possible that the amount of baryon asymmetry generated at the 
phase transition is 10^ — 10^ larger than what is observed. 

We solve the full quantum Boltzmann equations to extract the CP violating source 
terms which also provides enhancement compared to classical Boltzmann equations due to 
quantum memory effects . 

It is necessary to wash out some of the asymmetry; therefore the constraint v (Tc) /Tc ^ 1 
can be relaxed. We reevaluate the washout calculation using the effects of requiring some 
washout to occur and derive our limits for the light Higgs mass and the right handed stop 
based on electroweak baryogenesis. We are able to show that even in the one-loop approxi- 
mation for the temperature dependent effective Higgs potential the light Higgs mass upper 
limit can be pushed up to 115 GeV without invoking negative values of the right handed 
stop soft breaking mass parameter m^, which can lead to color-breaking global minima and 
scalar potential instability M. 
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II. BARYON ASYMMETRY CALCULATION 



In order to calculate the baryon asymmetry of the universe, one has to start with a 
self-consistent computation of the CP-violating sources resulting from particle interactions 
with the changing Higgs profile in the bubble wall. We follow the procedure of reference 
33| , which uses the closed-time-path formalism of finite-temperature field theory to derive 



quantum Boltzmann equations for Higgsinos and stops. Quantum memory effects resulting 
from correct treatment of particle propagation in the plasma strengthen the non-equilibrium 
character of the the particle scattering off the bubble wall and produce larger source terms. 

The dominant sources of CP-violation that are relevant to the electroweak baryogenesis 
scenario comes from the interactions of the Higgs fields with charginos and neutralinos. 
These interactions couple the Higgsino and gaugino components of charginos and neutralinos 
and involve potentially large CP-violating phases originating from the mixing. For the 
charginos, the CP-violating Lagrangian in the symmetric phase is 

C = -qH^HPlW - qHIWPlH + h.c. (1) 

The CP-violating phases ip^j, and (p2 are introduced by switching to the basis of mass eigen- 



states in the broken phase (see Appendix). Following the steps in to generate the 
CP-violating source term in the Boltzmann equation, we obtain for any point X inside the 
bubble wall 

Sc = 2gHm{M2fJ,y{X)P{X)I^ = sm{if^ + y^a), (2) 

where v"^ = vl + v^ and I3{X) = d(3{X)/dt characterizes the temporal variation of the Higgs 
profile tan/3(X) = t>2(X)/t>i(X) as the wall passes through point X. T^r is a tempera- 



ture dependent phase-space integral (the explicit form can be found in ||3^) and includes 
information about thermal behavior of the winos and Higgsinos in the high temperature 
plasma. 

The neutralino CP-violating interactions are 

C = -]^[HlH^PL{g2W^ - 9iB) + H^2{92W^ " 9iB)PlH2] + h.c. (3) 
Following the same steps as in the chargino case, we find the source term 

= gllm{M2fiy{X)P{X)I^ + gllu,{M^^i)v\X)i3{X)T^ 

= 7vFsin(¥'/. + </'2)+7Bsin((/?^ + ¥?i). (4) 

We can combine the source terms for the charginos and neutralinos to obtain the total 
source term for Higgsinos 

= 37vy sin(<^2 + Vt^) + 1b sm{ipi + if,,) (5) 



where 



7w = MM2\gyix)Pix)i^, (6) 

j^ = MM,\gy{X)${X)I^. (7) 
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Here we emphasize that the source dependends on the full physical (reparametrization 
invariant) combinations of the phases ipi + y?^ and (p2 + (fi/^ which factorize from the rest of 
the source term. In this sense our considerations are independent of the particular details 
going into the calculation of the Higgsino thermal production rate. 

The phase space integrals I^r and Xg exhibit strong resonant behavior leading to a 
maximum for mg rriyy [fn^j ~ m^) |^3[|. In terms of the soft breaking parameters the 



enhancement occurs when the gaugino masses Mi or M2 are close in value to the Higgsino 
mass parameter /i. For similar reasons it is easy to understand why the right handed stop 
contribution to the CP-violating source is always subdominant. The corresponding phase- 
space integral X^^ is typically far away from its maximum as m^^ ^ m^^ in our framework 
and the stop source can therefore be neglected. 

The baryon asymmetry can be directly related to the Higgsino source by solving a set 
of coupled diffusion equations for the Higgs and Higgsino, top and stop, and the first two 



generation quark and squark densities [0. We assume that strong sphaleron transitions 
and interactions induced by the top Yukawa coupling are very fast, allowing us to reduce 
the number of relevant equations to the diffusion equations for Higgses and Higgsinos and 
baryon number. They can be solved analytically, with the result 



(8) 



where D is an effective diffusion constant, T^^g is the weak sphaleron transition rate, and v^, 
is the velocity of the bubble wall. The coefficient A is determined by boundary conditions 
to be 

1 nOO 

A=jr^ duSn e-"+", (9) 
Jo 

where 



^ 2D ^ ^ 

and r is an effective decay constant. 

III. CALCULATION OF WASHOUT 

Once the phase transition has occurred, weak sphaleron transitions tend to erase any 



asymmetry that has been created. In the broken phase, the weak sphaleron rate is |28 



r^, ~ 2.8 X lO^T(^)^K(^ye^ 

where 



(12) 

aw 
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and 

B = B{^)c^l.87. (13) 

K is a constant in the range 10~^ < < 10^^. The differential equation satisfied by ub in 
the broken phase is 

= -CnfT^sriB, (14) 

where C is a number of 0{1) ||29| and we take C ~ 1 absorbing the uncertainty into the 
uncertainty of n and B. Solving (14) we get 

UBit) = nsitc) exp{- f dtUfV^s)- (15) 

The traditional bound is found by assuming that all of the washout happens at T = Tc, 
which amounts to saying that the integral in equation (3.5) is equal to the value of the 
integrand at T = T^. However, we want to evaluate the integral more carefully, and in the 
process derive a new bound. 

Since we do not have an analytic formula for v{T), we need to examine the integrand 
to see if we can make any simplifying assumptions. The sphaleron rate is proportional to 
{v(T)/Ty exp{—36v(T)/T), and we know that 1 < v(T)/T < oo. If we approximate that 
v{T) = v{Tc) throughout this period, then the washout rate is negligible when T ~ .75Tc ~ 
75 GeV. Since v{T) changes by about a factor of 2.5 in the temperature range from (0 - 100) 
GeV, treating v(T) as a constant is a good approximation; it is also conservative, because 
it will actually overestimate the washout slightly. 

We can make a change of variables to 

C = (16) 



by using the relation between time and temperature 

Mp 

The integral we now have to do is 



t (3 X 10-^)^. (17) 



/ = (3.4 X 10-')-=%-K r C'e-^dC 

Esp[Tc) JCc 

= (3.4xlO-)-^.e"^^xEl:C (18) 

Because Cc = '^'^Bv{Tc) / gTf. ~ 36, the Ql term will dominate all of the other terms. (Again, 
dropping the smaller terms is a conservative approximation.) Then 

/~(3.4xlO-«)^«:e-Hc 
~ (4.1 X lO^X.^e-"^^ (19) 
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Cc - 6 log Cc - log - 9 log 10 - log 4.1 + log(log /j, ~ 0- (21) 



Inserting this expression into equation (15), we now have an approximate expression for 
the observed ns/s: 

neiT ~ 0)/s = (^^(r = T,)/s) exp(-(4.1 x lO^X^e-^^) > 4 x (20) 
from which we can obtain 

nB/s{T c) 
4 X 10- 

By choosing a certain value for k and ub/ s{T^, this equation can be solved numerically. We 
will do this in the next section after we have presented our numerical result for nB/s{Tc). 
Once this is done, we can use the relation 

Cc ^ 36^ (22) 

^ c 

to find the lower bound on v{Tc)/Tc. 

This lower bound on v{T^/Tc can be translated into an upper bound on the light Higgs 
mass. The one-loop finite-temperature effective potential has a form 

V{^) = -\m\T)ip^ - T[EsMV^ + Fmssm{v.T)] + \\{T)ip\ (23) 

Z o 

For the MSSM with heavy decoupled stops, the potential becomes SM-like and one has [H 

v{Tc) _ 2EsM 
Tc ~ X 

where 



(24) 



2 M^ + M | 



= (25) 

Using ml = 2Af ^, equation becomes 



2 < 4:Esmv^ , . 

v{Tc)/Tc 

The other extreme is a light right-handed stop whose temperature dependent self energy, 
responsible for screening of the stop interactions in the plasma, is balanced by a negative 

soft squared mass term -mfj ~ U^^ = (|^| + ^g''^)T^ + i/i^ [1 + sin^ - ^)]T^. In such 
a case Fmssm is 

Fmssm {ip,T) = ip'EMssM (27) 

1-^ 2 

where At = At — ^/ tan/? is the stop left-right mixing parameter. Then 

vjTc) ^ 2{EsM + Emssm) 



(29) 



2 ^ 4:{EsAi + EmssmW 

When screening is present {mfj + U^^ > 0) the temperature dependent Higgs potential can 
be analyzed numerically. 
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IV. NUMERICAL RESULTS 



For the calculation of the baryon asymmetry, we have to adopt a model for the elec- 
troweak phase transition and evaluate v'^{X) and (3{X). These functions have been calcu- 
lated numerically using the two-loop finite-temperature effective potential in |3^. However, 
we have instead chosen to use the model of which lends itself more easily to our com- 
putation: 

v{X) = v{L^)[l - cos(^)][^^(X) - e{X - L^)] + v{L^)e{X - L^), (31) 
f3{X) = Af3[l - cos{^)][e{X) - e{X - L^)] + I3{X = 0) + A/50(X - L^). (32) 

We have tested that the answer only weakly depends on the model used for the phase 
transition, so the results should be of at least the right order of magnitude. We chose the 



value A/5 = .001 which is suggested by the results of ||36|, and a value v{L^) ~ 100 GeV. 

The width of the wall is chosen to be Zvji, ~ 25 /T [ PQI and the typical velocity of the bubble 
wall is expected to be ~ 0.1. For the averaged diffusion coefficient we use D = 0.8 GeV~^ 



and F = 1.7 GeV as in [R3|. All of these parameters enter the calculation of the baryon 



asymmetry, however, the scaling dependence of on these quantities is straightforward 
and does not interfere with CP-violation effects. 

Figure 1 shows the baryon asymmetry ub/s generated as a function of |/i| for several 
combinations of and (fi and if2 = 0, |Mi| = 140 GeV, IM2I = 250 GeV, tan/? ~ 3. The 
asymmetry hb/s can be as large as 10"'^ for ~ IM2I. 

Let us note in passing that the signs of ipi and most importantly of ip^ (as it appears 
in both chargino and neutralino matrices) determines the sign of n^. This fact is also 
demonstrated in Fig. 1 since ub is positive for siny^^ > with the exception of case e and 
fi ~ Ml where ub turns negative being dominated by the neutralino contribution to the 
Higgsino source. Experimental determination of the magnitude and sign of the soft phases 
is therefore essential if the feasibility of the supersymmetric electroweak baryogenesis is to 
be verified. Alternatively, the observed sign of could be used to determine the sign of 9?^ 
until it can be measured other ways. 

In order to calculate the lower bound on (c, we have to specify the the sphaleron parame- 
ters B and k. The non-perturbative scaling factor B comes from numerical minimization of 
the sphaleron energy and was evaluated for the MSSM in The usual range depending 



on the coupling strengths is 1.5 < -B(p-) < 2.7 with a typical median of 1.87. The value of 
K, is obtained as a functional determinant associated with fiuctuations about the sphaleron 
and was estimated to be k ~ 0.1 [0]. However, when the Higgs propagator uncertainty is 
absorbed the allowed range for kappa is 10~^ < k < 10~^. 

The required value of v{Tc)/Tc depends on the amount of dilution of the baryon number 
produced during the electroweak phase transition that can be allowed in order to explain 
the observed value of ub / s. Inclusion of both large CP- violating soft phases and quantum 
memory effects leads to a substantial enhancement in the value of the produced asymmetry 
nB/s{Tc). Consequently, the constraints on the strength of the first order phase transition 
can be softened as a function of nB(Tc)/nB{0). This decrease of v(Tc)/Tc is ilustrated in 
Fig. 2 for several values of k. The enhancement in nB{Tc) coming from quantum effects 
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FIG. 1. Plot of the baryon asymetry produced at the electroweak phase transition for 
\Mi\ = 140 GeV, IM2I = 250 GeV, tan/3 ~ 3 and f2 = 0. The five cases correspond to some 
typical cases, with a ip^ = ^, cpi = ^, h (p^ = ^, = |, c 99^ = ^, cpi = ^, d cp^ 
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10' 



"Pi 
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e ip^ = 5 X W (^1 = 5 X 10 



-3 



10' 



(about a factor of 10^ ||3^) for typical order of magnitude values of sin 99^ is demonstrated 
by the intervals with the position of the left (right) arrow corresponding to no (full) quan- 
tum enhancemnent respectively. The minimal value of sin required in order to generate 
nB/s{Tc) ~ nB/s{0) (negligible washout) without including the quantum effects in the 
sources is about 5 x 10~^ which is in agreement with the values obtained in Ref. [^. 

Since we are looking for an absolute lower bound on v{Tc)/Tc, we will take ns/siTc) to 
have its maximum value, 10"''. For k = 10~^, the numerical solution gives 



Cc > 39.6, 

while the solution for k, = 10~^ yields 

Cc> 31.3. 

This translates into a bound on v{Tc)/Tc of 

>.87, K = 10-\ 



(33) 
(34) 

(35) 
(36) 
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FIG. 2. Plot of required v{Tc)/Tc ratio as a function of the actual amount of baryon asymmetry 
produced at the electroweak phase transition. The intervals for different values of sin illustrate 
the enhancement resulting from quantum memory effects with the left (right) arrow corresponding 
to no (full) quantum enhancement. 



Using the bound for k = 10 ^, we can find the upper bound on the light Higgs mass. For 
heavy stops, the upper bound is 

m;,<51.7GeV. (37) 

For light stops with no thermal screening and neghgible mixing {Aj « mq) the upper 
bound is 



m/, < 138.1 GeV. (38) 

In Figure 3 we show a plot of the upper Higgs mass limit as a function of the right handed 
stop mass calculated in our framework. The stop mass range corresponds to variation of the 
right handed stop soft mass parameter roy from — (no thermal screening) to 11^^ (strong 
thermal screening). The wide band corresponds to variation of B and k within their full 
range and the narrow central band shows the set of curves resulting from taking B — 1.87 
and varying k in the full range. 

It is important to stress that the role of large CP-violating phases is crucial in this 
context. If the baryon asymmetry is overproduced during the electroweak phase transition 
due to 0{1) values of ip^ and/or tpi the washout conditions (35) and (36) are less stringent 
then if the phases are constrained to be 10~^. As a result, the upper bound on the Higgs 



10 



180.0 
170.0 
160.0 
150.0 
140.0 
130.0 
120.0 
110.0 
100.0 
90.0 



80.0 



130.0 



1.5<B<2.7 

10"^<K<10'' 




150.0 



170.0 



190.0 



210.0 



FIG. 3. Dependence of the upper Higgs mass bound on the light right handed stop mass. The 
width of the bands results from uncertainties in the weak sphaleron parameters B and k which 
are varied in the full range. The central band corresponds to B=1.87 while k is varied in the full 
range. 



mass can be as much as 15 GeV higher compared to the situation where nB/s{Tc) ~ nB/s{0) 
and no washout is allowed. 

It is obvious from our results that inclusion of large CP- violating phases in the calculation 
of the baryon asymmetry relaxes the stringent constraints on the strength of the first order 
transition and consequently the light Higgs masses can be pushed towards larger values. 
Even for stop masses m^^ ^170 GeV resulting from positive values of the soft breaking 
parameter mfj the upper bound on the Higgs mass can be as high as 115 GeV. 

Our considerations are based on the one loop temperature dependent effective Higgs 
potential evaluation. Two loop corrections are known to significantly enhance the strength 



of the first order phase transition [|I^ and further relax the upper bound on the Higgs mass. 
In this respect our results represent a conservative estimate of the upper Higgs mass limit 
and it is likely to be moved upward when two loop corrections are included. 



V. SUMMARY AND CONCLUSIONS 

The upper bound on the Higgs mass which still allows electroweak baryogenesis is an 
important issue which is becoming relevant as the Higgs mass experimental lower limits 
are approaching 100 GeV. The baryon asymmetry produced during the electroweak phase 
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transition can succesfully explain the observed value provided there are additional degrees 
of freedom contributing to the finite temperature Higgs potential. Supersymmetric models 
with a light right handed stop are a very good candidate for a theory that can provide these 
degrees of freedom. Also, the CP-violating phases appearing in the soft breaking terms of 
the Lagrangian can supplement (or entirely replace) the effects of a CP-violating phase in 
the CKM matrix. 

Previously it has been thought that based on the electron and neutron EDM experimental 
limits the supersymmetric CP-violating phases have to be small and consequently there is 
no room for washout of the produced baryon asymmetry. This translates into stringent 
constraints on the Higgs and right handed stop masses, often leading to problems with color 
breaking minima and potential stability. 

We have shown that once there is a possibility of cancellations among individual con- 
tributions to the EDMs and the CP- violating phases are allowed to be 0{1), the produced 
baryon asymmetry exceeds the observed value and the baryon density can be allowed to be 
diluted by a factor of 10"^ or more. In such a case the upper Higgs mass bound is increased 
and depending on the right handed squark mass it can go beyond 115-120 GeV while at the 
same time the scalar potential is stable. 

This scenario opens new possibile implications for supersymmetric phenomenology. As 
pointed out in it is important to independently measure the CP-violating phases to 



correctly interpret experimental observables. If the phases are small the window for super- 
symmetric electroweak baryogenesis will get increasingly smaller as the LEP and Tevatron 
experiments will be pushing up the lower Higgs mass limit. It is natural to expect that in 
the case of small phases the experimentally determined Higgs mass should not be too far 
above 100 GeV if electroweak baryogenesis is expected to work. On the other hand, if the 
Higgs is not discovered at LEP nor at the Tevatron, one can expect that the CP-violating 
phases of the MSSM are of 0{1) if baryogenesis occurs at the electroweak phase transition, 
and they should be measurable if the superpartners are discovered. Of course, finding a 
light Higgs boson at LEP or Fermilab is consistent with having large supersymmetric soft 
phases. 
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VII. APPENDIX 

Above the electroweak scale, the chargino mass matrix is 

/ iMsle^^^ \ ,,,, 
^-^[ Me- J- 

This matrix is made real and diagonal by two complex matrices. 
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where we can take 



U 



1 
1 



e^'^'^ 



If we switch to a basis of two-component Weyl spinors, 

PlW = PlV0j, PrW = PrUji^Ijj, 

PlH = PlV0j, PrH = PrUj2^I^j, 
where ijjj are two-component spinors. In terms of ipj, the interaction becomes 
£ = -gHle-''P^tl]2PLA - gH^e-^'^'^ilJiPLih + h.c. 
Above the electroweak scale, the neutralino mass matrix is 



Mn — 



/|Mi|p>i 

|M2|e''^2 



V -|^|e*^^* 



\ 


J 



It is made real and diagonal by the complex matrix N, 

Mt""' = N^MnN, 

where we can take 





TV 












V 



4.e%^ ^e^ 



V2 

Switching to the two-component spinor basis, we have 



(A2) 
(A3) 

(A4) 
(A5) 

(A6) 
(A7) 
(A8) 



(A9) 



PlHi 




PrHi 


= PrNj^^i^i 


(AlO) 


PlH2 


= PlN^^iPI 


PrH2 


= PrNj4^PI 


(All) 


PlWs 


= PlN*Ji 


PRm 


= PRNJ2^P^, 


(A12) 


PlB 




PrB-- 


-- PrNji^I 


(A13) 



resulting in the interaction 



(A14) 
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